First principles study of the crystal structure and dehydrogenation pathways of Li4BN3Hio 
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Using density functional theory we examine tlie crystal structure and tlie finite-temperature thermodynamics 
of formation and dehydrogenation for the new quaternary hydride Li4BN3Hio. Two recent studies based on 
X-ray and neutron diffraction have reported three bcc crystal structures for this phase. While these structures 
possess identical space groups and similar lattice constants, internal coordinate differences result in bond length 
discrepancies as large as 0.2 A. Geometry optimization calculations on the experimental structures reveal that 
the apparent discrepancies are an artifact of X-ray interactions with strong bond polarization; the relaxed struc- 
tures are essentially identical. Regarding reaction energetics, the present calculations predict that the formation 
reaction 3 LiNH2 -I- LiBH4 Li4BN3Hio is exothermic with enthalpy A_ff'^='"""^ = -11.8 kJ/(mol f.u.), con- 
sistent with reports of spontaneous Li4BN3Hio formation in the literature. Calorimetry experiments have been 
reported for the dehydrogenation reaction, but have proven difficult to interpret. To help clarify the thermody- 
namics we evaluate the free energies of seventeen candidate dehydrogenation pathways over the temperature 
range T — 0-1000 K. At temperatures where H2-release has been experimentally observed (T « 520-630 K), 
the favored dehydrogenation reaction is Li4BN3Hio Li3BN2 -I- LiNHa -I- 4H2, which is weakly endothermic 
[AH^^^^^^ — 12.8 kJ/(mol H2)]. The small calculated AH is consistent with the unsuccessful attempts at re- 
hydriding reported in the literature, and suggests that the moderately high temperatures needed for H-desorption 
result from slow kinetics. 



PACS numbers: 64.70.Hz, 65.40.-b, 71.15.Nc, 81.05.Zx 

I. INTRODUCTION 

Recent efforts to improve the efficiency and reduce the en- 
vironmental impact of automobile transport have focused on 
hydrogen-based fuel cells (FC) and internal combustion en- 
gines (H2ICE) as possible replacements for current technolo- 
gies powered by fossil-fuels.' A significant obstacle to realiz- 
ing this transition is the on-board storage of hydrogen at high 
gravimetric and volumetric densities.^ To achieve the storage 
densities necessary for mobile applications, novel means for 
H-storage are necessary, and a set of targets have been estab- 
lished to guide the search for new storage systems.^ At present 
no known storage material or mechanism meets these targets. 

One promising avenue for efficient storage of hydrogen is 
via solid state storage, such as in the form of complex- or 
metal-hydrides. Solid storage has the advantage of provid- 
ing volumetric densities beyond what can be achieved with 
compressed H2 gas or cryogenic liquid storage. However, all 
known hydrides suffer from one or more of the following lim- 
itations: low gravimetric densities, high H-desorption temper- 
atures, or an inability to be easily re-hydrided. The search 
for new hydrides that overcome these limitations has attracted 
intense interest during the past 5 years. 

Towards these ends, recent experiments by Pinkerton and 
co-workers"^'^'^ and Aoki et alr^ on the H-storage properties 
of the quaternary Li-B-N-H system are noteworthy. By mix- 
ing lithium borohydride (LiBH4) and lithium amide (LiNH2), 
both groups have reported the formation of a new hydride 
phase, which when heated above '^520 K released approx- 
imately 10 wt.% hydrogen. While the reported desorption 
temperature is somewhat too high, and reversibility has not 
yet been demonstrated, improvement in these areas may be 



possible (and has been partially demonstrated'') through the 
addition of catalytic dopants^-^ or via novel synthetic routes,— 
suggesting that further study of this system is warranted. 

The composition of the new Li-B-N-H phase was prelimi- 
narily identified"* as "Li3BN2H8," but subsequent experiments 
based on single crystal X-ray diffraction," and, synchrotron 
X-ray and neutron powder diffraction, identified its true stoi- 
chiometry as Li4BN3Hio. Combined, these diffraction experi- 
ments identified three similar crystal structures, all sharing the 
bcc space group /2i3 with lattice constant a = 10.66-10.68 A, 
but with somewhat different internal coordinates. The dis- 
crepancy in atomic positions results in N-H bond lengths that 
differ between structures by as much as 0.2 A, with smaller 
N-H lengths reported for the X-ray structures. It has been 
suggested that the discrepancy may be an artifact of X-ray in- 
teractions with strongly polarized electron clouds along the 
N-H bond direction»ii 

A careful characterization of a hydride's crystal struc- 
ture is desirable because it enables an independent, ab ini- 
tio assessment of the thermodynamics of hydrogen desorp- 
tion/absorption reactions. Such an assessment is of value 
because it can clarify the thermodynamics when calorimetry 
measurements yield ambiguous results, such as in situations 
where multiple reactions occur simultaneously (see below). A 
key thermodynamic property which determines the suitabil- 
ity of a hydride for H-storage applications is the strength of 
the hydrogen-host bond. The bond strength is quantified via 
the change in enthalpy {AH) occurring during hydrogen up- 
take/release. In mobile FC applications, for example, it is de- 
sirable that the desorption reaction be ent/of/ierm/c {AH > 0) 
with AH w 20-50 kJ/(mol H2). Neglecting kinetics, an en- 
thalpy in this range would allow for H-desorption at temper- 
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atures compatible with the waste heat from a FC, and permit 
on-board recharging under reasonable H2 pressure. 



Unfortunately, calorimetric measurements of H-desorption 
in Li4BN3Hio have been difficult to interpret. (See Refs.|6jand 

and the supplementary information accompanying Ref. 
Because Li4BN3Hio melts at ^460 K, approximately 60 de- 
grees below the onset of H-desorption, hydrogen gas evolves 
from the liquid stateii Concurrent with H2 release is the for- 
mation of a solid reaction product, generally involving one 
or more polymorphs of Li3BN2 along with other unidentified 
phases, ^'^'^ and desorption of 2-3 mol % ammonia."* Measure- 
ments based on differential scanning calorimetry (DSC)i^ and 
differential thermal analysis (DTA)** suggest that the net heat 
flow during dehydriding is exothermic (AiJ < 0). However, 
it is unclear what fraction of the total thermal profile arises 
from the exothermic latent heat of Li3BN2 solidification,'^ and 
furthermore how H2 and NH3 release impacts the calorime- 
ter response.'* Consequently, it has not been possible to un- 
ambiguously assess the endo- or exothermic nature of H- 
desorption alone. Neglecting kinetic effects, the fact that it 
has proven very difficult"* to re-hydride Li4BN3Hi() suggests 
desorption is weakly endothermic, or exothermic. Regard- 
ing Li4BN3Hio formation, several studies^'^ '^ have reported 
the spontaneous formation of Li4BN3Hi() after mixing LiBH4 
with LiNH2, even at temperatures as low as room temperature, 
suggesting that the formation reaction is exothermicj^ 



In light of the discrepancy noted above regarding the crys- 
tal structure of Li4BN3Hi(), and the additional ambiguity sur- 
rounding its dehydrogenation thermodynamics, we employ 
density functional theory (DFT)*^'*"* calculations in an attempt 
to clarify these issues. First, we determine the ground state 
crystal structure by performing separate geometry optimiza- 
tion calculations on each of the experimentally proposed crys- 
tal structures. We find that the relaxed structures are essen- 
tially identical, and have N-H bond lengths that are consistent 
with the structure based on neutron diffraction,'^ thus con- 
firming the conjecture" that the anomalously short N-H bond 
lengths observed in Ref. 11 are an artifact of X-ray measure- 
ments. Second, we evaluate the finite-temperature reaction 
enthalpies and free energies for the formation and dehydro- 
genation of Li4BN3Hio. Li4BN3Hi() formation (with respect 
to LiNH2 and LiBH4) is found to be exothermic, in agree- 
ment with experimental reports of spontaneous Li4BN3Hio 
formation in the literature. '*'^'*^ For H2-desorption we ex- 
plore the thermodynamics of seventeen candidate reactions 
over the temperature range T — 0-1000 K, as there appears 
to be some uncertainty in determining the reaction products 
experimentally.'*'^ At temperatures where desorption has been 
reported"* {T k, 520-630 K), the favored reaction is predicted 
to be Li4BN3Hio —> Li3BN2 + LiNH2 + 4 H2, which is weakly 
endothermic. Moreover, the calculated free energies suggest 
Li4BN3Hio is a metastable phase that should decompose via 
one of three temperature-dependent pathways. 



II. METHODOLOGY 

First-principles calculations were performed using a 
planewave method based on the PW91 generalized gradient 
approximation'^ to density functional theory'-'-''* (VASP).'- 
The core-valence electron interaction was treated using 
Blochl's projector augmented wave (PAW) method,'^ '^ and k- 
point sampling was performed on a dense Monkhorst-Pack^^ 
grid with an energy convergence of better than 1 meV per su- 
percell. (In the case of Li4BN3Hi(), a 2 x 2 x 2 grid was 
sufficient to acheive this level of precision.) Electronic occu- 
pancies were determined via a gaussian smearing algorithm 
with a 0.1 eV smearing width. For high-precision calcula- 
tion of static zero Kelvin (omitting zero point vibrational ef- 
fects) electronic energies and crystal structures we used the 
so-called "hard" VASP PAW potentials^' with a planewave cut- 
off of 875 eV, and a geometry relaxation tolerance of better 
than 0.02 eV/A. Internal atomic positions and external cell 
shape/volume were optimized simultaneously. 

Finite-temperature thermodynamics were evaluated within 
the harmonic approximation.^"-^' Vibrational frequencies {oji) 
were extracted by diagonalizing a dynamical matrix whose el- 
ements were determined via the so-called direct method:^^ the 
forces generated by series of symmetry inequivalent atomic 
displacements about the equilibrium geometry (±0.02 A, 
±0.04 A) were fit to cubic splines in order to extract the 
force constants. Because large supercells are desirable to min- 
imize finite-size effects, and with the exception of the molec- 
ular species (H2, N2, NH3) where the number of atoms per 
supercell is small, the dynamical matrix calculations on the 
solid-state phases were preformed using a softer set of PAW 
potentials'^ and planewave cutoff energies ranging from 400- 
500 eV.'2 A comparison of calculated structural parameters 
and formation energies with existing experimental and theo- 
retical data is presented in the following section. 

Once the normal-mode frequencies have been determined, 
finite-temperature energetics can be obtained by enthalpic 
(i/vib) and entropic (S'vib) additions to the static electronic en- 
ergies. Within the harmonic approximation these contribu- 
tions are given byi^ 



i?vib(T) 



explTT^j - 1 
Kb -I 



1 -1 



(1) 



■In 



1 — exp(' 



(2) 



where the sums run over vibrational frequencies (3N — 3 fre- 
quencies for solids and non-linear molecules, 3N — 5 fre- 
quencies for linear molecules), fc^ is the Boltzmann factor, 
and T is the absolute temperature. For the linear (non-linear) 
molecules an additional ^ksT (iksT) term is added to Eq.[T] 
to account for translational, rotational, and pV degrees of free- 
dom. The zero point energy (ZPE) can be recovered from 
Eq.[T]in the limit H^ib{T = 0). The enthalpy and free energy 
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of a phase can therefore be expressed as: 

H{T) = E + H,,^{T) (3) 
G{T) = H{T)-S(T)T, (4) 

where E is the static electronic energy of the crystal/molecule 
in its ground-state geometry and S represents either the stan- 
dard tabulated-^ entropy of a given molecular species in the 
gas phase at p = 1 bar [S^=^^^^ = 130.858 (H2), 191.789 
(N2), and 192.995 J K'^ mol"! (NH3)], or the vibrational en- 
tropy S'vib of a soUd state phase. 

III. CRYSTAL STRUCTURE 

The crystal structure of Li4BN3Hio has been studied by two 
groups." '^ Filinchukand co-workers" used single-crystal X- 
ray diffraction to identify the structures of the two largest crys- 
tal domains obtained after re melting a 2:1 mixture of LiNH2 
and LiBH4. Both domains exhibited a bcc structure (space 
group /2i3) with lattice constants of 10.67-10.68 A. Chater et 
al— used a combination of high-resolution synchrotron X-ray 
and neutron diffraction to examine powder samples prepared 
from a wide range of LiNH2:LiBH4 compositions. They de- 
termined that although the most likely stoichiometric compo- 
sition was 3:1, the Li4BN3Hio structure was able to accom- 
modate a wide range of stoichiometrics. Their best-fit crystal 
structure had external cell parameters similar to those reported 
in Ref . , 1 1; bcc (space group /2i3) with a — 10.66 A. Despite 
the good agreement in external geometries, the structures ob- 
tained by these two studies differ markedly in their internal 
coordinates. Most notably, the single-crystal X-ray diffrac- 
tion study" found anomalously short N-H bond lengths of 
0.83-0.86 A, while neutron scattering'^ gave lengths of 0.98- 
1.04 A. 

The external cell parameters obtained upon relaxing each 
of the three Li4BN3Hio experimental structures (as well as the 
structures of the other phases used in our subsequent discus- 
sion of reaction thermodynamics) are summarized in Table J] 
The resulting structures are bcc and share the same lattice con- 
stant, a — 10.60 A, in good agreement with both diffraction 
studies. 

For each of the three relaxed structures. Table com- 
pares the calculated internal atomic coordinates with the cor- 
responding experimental values. We note first that the aver- 
age absolute deviation (5) between theory and experiment is 
smallest for the neutron structure,'^ 5 — 3.6 x 10^"^, whereas 
for the X-ray sti-ucture" 5 = 4.4 x 10"^ and 5 = 4.2 x 10"^ 
for domains 1 and 2, respectively. The largest discrepancy 
between theory and experiment is in the positions of the N- 
bonded hydrogen atoms, HI & H2, in the X-ray structure. 
(Note that the labeling convention for hydrogens is different 
in the neutron structure: there HI & H2 refer to B-bonded hy- 
drogens.) Secondly, while there appear to be large differences 
in the internal coordinates measured by Refs. ITU and iT2l (Ta- 
ble HHi. after relaxation these three structures may be mapped 
onto one another via a series of rigid body rotations and trans- 
lations. We conclude that all three experimental structures re- 
lax to essentially the same structure. 



TABLE \: Calculated structural parameters compared with experi- 
mental data. Bond lengths (d) and crystal lattice constants (a, b, c) 
are given in A, bond angles (Z) and angles between lattice vectors 
(q, P) are listed in degrees. Space group information for the solid 
phases is listed in parentheses. 



System 


Parameter 


Calculated 


Experiment 


H2 


d(H-H) 


0.749 


0.741^ 


N, 


d(N-N) 


1.101 


1.098ii 


NH3 


d(N-H) 


1.021 


1.0122i 




Z(H-N-H) 


106.4 


106.7^4 


a-B (_R3m) 


a 


5.05 


5.06^^ 




a 


58.1 


58.1^^ 


BN (F43m) 


a 


3.62 


3.632^ 


Li (/m3m) 


a 


3.44 


3.51.2i 


LiH (FmSrn) 


a 


4.01 


4.01^ 


LisN (P3ml) 


a 


3.64 


3.61° 




c 


3.88 


3.85° 


LiNH2 (/4) 


a 


5.02 


5.04^^ 




c 


10.28 


10.28-' 


Li2NH (Pnma) 


a 


7.75 


7.73' 




b 


3.61 


3.60* 




c 


4.88 


4.87* 


a-Li3BN2 {P4:2/mnm) 


a 


4.67 


4.64'°=^ 




c 


5.22 


5.26'°^ 


I3-U3BN2 (P2i/c) 


a 


5.15 


5.I5S 




b 


7.08 


7.08^ 




c 


6.77 


6.79S 




P 


112.7 


113.0^ 


Li3BN2 (Mi/amd) 


a 


6.63 


6.6(^ 




c 


10.35 


10.3^ 


LiBH4 (Pnma) 


a 


7.25 


7.18^ 




b 


4.38 


4.44^ 




c 


6.63 


6.80^ 


Li4BN3Hio (/2i3) 


a 


10.60 


10.67-10.68- 

io.6ei^ 



"Theoretically predicted structure from Ref. 35. In agreement with Ref. [33. 
we find the room-temperature experimental P6 / mmm structure to be unsta- 
ble at low temperature, with a weak soft mode of 66i cm~^. 

''Theoretically predicted structure from Ref. [33. 



A comparison of calculated bond lengths, bond angles, and 
intermolecular distances is given in Table [III] Overall, the 
relaxed structures agree best with the /nframolecular lengths 
(i.e., within an NH2 or BH4 fragment) determined using neu- 
tron diffraction; conversely, the /nfermolecular distances from 
X-ray are closest to those in the relaxed theoretical structures. 
As for the intramolecular structure, as noted above, the major 
discrepancy between the experimental structures lies in the 
two N-H bond lengths, (i(N-H), in the NH2 fragment. The 
calculated length of 1.027 A is in better agreement with the 
longer bond lengths predicted by neutron diffraction of 0.983 
and 1.042 A. Similarly, the X-ray data also underestimates 
B-H bond lengths in the BH4 units relative to our DFT calcu- 
lations and neutron data. Calculated bond angles agree best 
with those from the neutron structure. 
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TABLE II: Calculated relaxed internal atomic positions of Li4BN3Hio compared with experimental measurements from Refs. fTTI I T^ . 







Calculated 






Experiment 






Atom 


X 


y 


z 


X 


y 


z 




N 


0.117 


0.362 


0.402 


0.115 


0.359 


0.405 




B 


0.113 


0.113 


0.113 


0.114 


0.114 


0.114 




Lil 


0.281 


0.000 


0.250 


0.287 


0.000 


0.250 




Li2 


0.519 


0.000 


0.250 


0.524 


0.000 


0.250 


Domain 1, 


Li3 


0.484 


0.484 


0.484 


0.483 


0.483 


0.483 




HI 


0.195 


0.306 


0.393 


0.175 


0.310 


0.397 


Ref.JJ. 


H2 


0.118 


0.415 


0.321 


0.115 


0.405 


0.339 




H3 


0.003 


0.118 


0.147 


0.012 


0.119 


0.146 




H4 


0.180 


0.180 


0.180 


0.172 


0.172 


0.172 




N 


0.117 


0.361 


0.402 


0.115 


0.359 


0.405 




B 


0.113 


0.113 


0.113 


0.114 


0.114 


0.114 




Lil 


0.281 


0.000 


0.250 


0.286 


0.000 


0.250 




Li2 


0.519 


0.000 


0.250 


0.523 


0.000 


0.250 


Domain 2, 


Li3 


0.484 


0.484 


0.484 


0.484 


0.484 


0.484 




HI 


0.195 


0.305 


0.394 


0.176 


0.310 


0.399 


Ref.JJ. 


H2 


0.118 


0.415 


0.321 


0.114 


0.402 


0.338 




H3 


0.003 


0.117 


0.148 


0.015 


0.116 


0.149 




H4 


0.180 


0.180 


0.180 


0.174 


0.174 


0.174 




N 


-0.152 


0.112 


0.367 


-0.156 


0.110 


0.367 




B 


0.137 


0.137 


0.137 


0.135 


0.135 


0.135 




Lil 


0.250 


-0.030 


0.000 


0.250 


-0.043 


0.000 




Li2 


-0.268 


0.000 


0.250 


-0.262 


0.000 


0.250 




Li3 


0.266 


0.266 


0.266 


0.272 


0.272 


0.272 


Ref.J2 


HI 


0.071 


0.071 


0.071 


0.072 


0.072 


0.072 


H2 


0.103 


0.133 


0.247 


0.097 


0.136 


0.248 




H3 


-0.143 


0.055 


0.445 


-0.153 


0.053 


0.439 




H4 


-0.071 


0.165 


0.368 


-0.077 


0.167 


0.377 





TABLE III: Comparison of calculated atomic distances and angles with values determined from single-crystal X-ray diffraction (Ref.fTlh and 
synchrotron X-ray and neutron powder diffraction (Ref. 12). Lengths (d) are given in A, angles (/) are in degrees. 



Ref. 11, domain 1 Ref. 11, domain 2 Ref. 12 



Parameter 


Calculated 


Experiment 


Calculated 


Experiment 


Calculated 


Experiment 








Intramolecular Values 








d(N-H) (0 


1.027 


0.83(2) 


1.027 


0.84(2) 


1.027 


0.983 


d(N-H) (ii) 


1.027 


0.859 


1.026 


0.852 


1.027 


1.042 


d(B-H) (0 


1.223 


1.08(3) 


1.225 


1.11(3) 


1.224 


1.169 


d(B-H) (n)x3 


1.225 


1.137 


1.225 


1.121 


1.223 


1.270 


/(H-N-H) 


103.5 


106.2 


103.5 


105.9 


103.5 


104.9 


/(H-B-H)x3 (0 


108.0 


108.9 


107.9 


109.2 


107.9 


107.6 


/(H-B-H)x3(h) 


110.9 


110.1 


111.0 

Intermolecular Values 


109.8 


111.0 


111.3 


d(N-Li) (i) 


2.071 


2.069 


2.069 


2.068 


2.070 


2.033 


d(N-Li) (ii) 


2.116 


2.115 


2.111 


2.113 


2.110 


2.053 


d(N-Li) (Hi) 


2.138 


2.117 


2.138 


2.119 


2.142 


2.094 


d(N-Li) (iv) 


2.164 


2.157 


2.164 


2.157 


2.169 


2.209 


d(B-Li) (i) 


2.371 


2.410 


2.375 


2.406 


2.368 


2.527 


d(B-Li) (i0x3 


2.595 


2.651 


2.592 


2.642 


2.589 


2.681 



IV. REACTION ENERGETICS 



In order to examine the finite-temperature thermodynamics 
of reactions involving Li4BN3Hio, it was first necessary to cal- 
culate the ground-state structures of the phases that participate 
in those reactions. A summary of the compounds used, their 



crystal structures, and a comparison of the calculated struc- 
tures to experimental data is presented in Table U In general 
the agreement between experiment and theory is very good. 

Some of the compounds listed in Table U are known 
to have several polymorphs. For example, at least three 
polymorphs have been reported for Li3BN2," ^''^^'^^'^^ includ- 
ing low-temperature tetragonal (a),"" high-temperature mon- 
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Temperature (K) 



FIG. 1: Calculated Gibbs free energies (in kj/mol f.u.) of mono- 
clinic (/?) and body-centered tetragonal (BCT) LIsBNt as a function 
of temperature. The energy zero is set to the static K energy of 
P-U3BN2. 



oclinic and high-pressure body-centered tetragonal 

(bct)S2iIL^ phases. All three of these phases have been ob- 
served as dehydrogenation products forLi4BN3Hio,'*'** and we 
have performed calculations on each phase to determine its 
relative stability as a function of temperature. The /3-Li3BN2 
phase was found to have the lowest zero-Kelvin static en- 
ergy, « 2 kJ/(mol f.u.) lower than that of either the bet or 
a phases. The latter two phases are degenerate in energy to 
within 0. 1 kJ/(mol f.u.). Our zero-Kelvin energetics and struc- 
tures agree well with recent calculations reported by Pinkerton 
and Herbst.^*^ 

Vibrational contributions to the energies of each Li3BN2 
polymorph were evaluated within the harmonic approxi- 
mation. For Q;-Li3BN2 a careful evaluation of the vibra- 
tional spectra yielded doubly-degenerate imaginary modes at 
81« cm^^, indicating that this phase is unstable at low temper- 
atures. A search for alternative low-energy a-Li3BN2 struc- 
tures was performed using molecular dynamics (md) on an en- 
larged (2x2x2) a-Li3BN2 supercell at T = 423 K for 10 ps. 
At 1 ps intervals the cuiTent md configuration was stored, re- 
laxed, and the symmetry of the resulting optimized structure 
was determined. Four unique crystal structures were identi- 
fied, with space groups (numbers): P2i2i2i (19), Pccn (56), 
Pmmn (59), and Pnraa (62). These structures were then 
relaxed within their symmetry constraints, and were found 
to have energies ^1 kJ/(mol f.u.) lower than the original 
a-Li3BN2 structure [which is still 1 kJ/(mol f.u.) higher in 
energy than /?-Li3BN2]. In light of these results we suggest 
that the crystal structure of a-Li3BN2 be experimentally re- 
assessed, and we exclude this phase from our evaluation of 
Li4BN3Hio dehydrogenation reactions. Further information 
regarding our proposed structure of a-Li3BN2 can be found 
elsewhere.''^ 

Unlike a-Li3BN2, no imaginary modes were observed in 
the vibrational spectra of (3 or bet Li3BN2. Using Eqs. [T]|2] 
Fig. [T] plots the Gibbs free energies (Eq. |4]i of these phases 
as a function of temperature. Due to their similar vibrational 



TABLE IV: Calculated vibrational contributions (at T = iQQK) to 
the free energies of phases used in this study. ZPE refers to the zero 
point energy, given by H^ibiT = 0) (Eq. [T); iSvib = ^^vib — ZPE 
(omitting the ksT molecular terms); Svib is the vibrational entropy, 
Eq.|2] Units are kj moP^ for ZPE and E^it,, and J moP^ K^^ for 



5'vib- 




ZPE 


7PF other 


TP T=300K 


Q T=300K 
'^vib 


ni 


25.7 


25.7^, 26.1—, 


« 


~ 












Nt 




14442 






i>iri3 


/ .0 


RQ 1^ 


n 1 


n A 




^ 
j.y 


J. 5*" 


4. J 


17 7 
Z / . / 


B 


12.5 


I2.3ii 


I.I 


5.3 


BN 


29.9 


30.9^ 






LiH 


21.7 


21.4,ii21.8^ 


3.5 


17.7 


LiaN 


28.0 


28.62^ 






LiNH2 


69.0 


69.1," 69.5, - 


8.9 


51.6 






69.8i^ 






Li2NH 


46.7 


46.7'*^47.1^ 


9.6 


52.9 


LiBH4 


107.1 


103.0,"' 106.5,- 


10.8 


63.6 






108. li^ 






Li3BN2 


52.2 




15.3 


83.8 


Li4BN3Hio 


314.4 




37.0 


213.4 



contributions (Table IIVI ). the free energies of the /? and bet 
phases exhibit nearly identical temperature dependence. Con- 
sequently, the small ^^2 kJ/(mol f.u.) difference in static ener- 
gies noted above is sufficient to favor the monoclinic f3 phase 
as the ground state structure, and the bct-/5 energy difference 
remains roughly constant for T < 1000 K. Based on the vi- 
brational instability observed in a-Li3BN2, and the higher free 
energy of bet Li3BN2 (relative to /?-Li3BN2), our subsequent 
thermodynamic analyses are performed assuming Li3BN2 re- 
sides in the the low-energy (3 structure. 

The zero-point, enthalpic, and entropic contributions to the 
free energy of the various phases are summarized in Table IIVI 
As a further test of our computational methodology, there we 
also compare our calculated ZPE to other values reported in 
the literature, and find good agreement. With the exception of 
Li, vibrational effects are found to contribute substantially to 
the free energies of these systems. For example, in Li4BN3Hio 
the ZPE exceeds 300 kJ/mol. 

A comparison of formation energies for the compounds 
used in this study (relative to the elements in their standard 
states) is presented in Table [V] Three reports^*^'"^'-'*^ based on 
DFT calculations have also recently evaluated formation en- 
ergies for some of these compounds, and the present calcu- 
lations are in good agreement with their findings. Further- 
more, the agreement with experimental enthalpies is reason- 
able (generally on the order of 10%), which is representa- 
tive of the accuracy obtained via density functional methods 
for reactions involving molecular species. As expected, the 
largest percentage discrepancy is for the ammonia molecule.— 
Two general trends in the data are evident: (/) With the excep- 
tion of ammonia, the DFT enthalpies (AiJ^^^*^'"^) are more 
positive than the experimental values, and (//) the inclusion of 
dynamical contributions (AE Ai/^^-'""^) results in more 
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TABLE V: Comparison of calculated formation energies (with re- 
spect to the elements, in kj/mol f.u.) with existing experimental and 
theoretical data. AE refers to differences in static K energies, for- 
mation enthalpies (AH) and free energies (AG) are based on Eqns.[3] 
and|4l respectively. Experimental data is from Refs.l23land[^ theo- 
retical data is from a set of recent DFT calculations (Refs. [33pl1l42l) . 



System 


AE 


^^T=300K 


^^T=300K 


NH3 


-98.0 


-63.2 


-33.6 






-45.9" 


-16.32^ 


BN 


-233.8 


-227.4 


-198.9 






-250.9^ 


-224.9" 


LiH 


-83.6 


-83.8 


-61.2 




-81^1- 


-90.62^ 


-68.3^ 




-83.9^^ 


-84.6^^ 




LigN 


-148.7 


-145.9 


-109.8 






-164.62^ 


-128.4" 


LiNHo 


-196.6 


-172.6 


-111.7 




-196.^ 


-176^ 








-173.li^ 




Li^NH 


-196.0 


-184.5 


-135.3 




-194.C- 


-222^^ 








-184.li^ 




LiBH4 


-205.9 


-178.6 


-109.2 




-194''!. 


-190.5^^ 


-124.32^ 


Li3BN2 


-497.4 


-490.6 


-433.6 




-49733 







positive formation energies. 

Using the energetic contributions from Table |IV| in Ta- 
ble [Vl] we evaluate the thermodynamics for two Li4BN3Hi() 
formation reactions. Reaction (a) considers formation from 
the elemental phases in their standard states, and reac- 
tion (b) corresponds to the synthesis route employed in the 
literature^'^ '- involving a mixture of LiNH2 and LiBH4. (We 
consider only the stoichiometric'^ 3:1 ratio of LiNH2 to 
LiBH4.^-) The energetics of each reaction are are decomposed 
into a sequence having increasingly more physical contribu- 
tions to the reaction thermodynamics, AE AH^^'^^ 
^^T=30()K _^ ^(^T=30()K^ which allows us to gauge the im- 
portance of these contributions in comparison to the common 
practice of evaluating only static zero-Kelvin energetics. [AE 
refers to differences in the static zero-Kelvin energies (E in 
Eq.O, and AH'^=°^ adds ZPE contributions to AE.] Both 
formation reactions in Table [Vl] are predicted to be exother- 
mic (AH < 0), in agreement with reports of spontaneous 
Li4BN3Hio formation from the literature. "''^ '^ As expected, 
formation from the elements vs. from the hydride mixture is a 
substantially more exothermic process. 

It is noteworthy that dynamical contributions to the reaction 
energies are significant in reaction (a) {AE ^ AH), but ap- 
pear to play a minor role in reaction (b) {AE w AH). By ex- 
amining the values presented in TablellVlfor the phases partic- 
ipating in reaction (b) (Li4BN3Hio, LiNH2, & LiBH4) we see 
that these terms are sizable, but that they largely cancel out. 
This cancellation effect can be understood by the similarities 
in internal bonding shared by these three phases: Li4BN3Hi(), 
with its Li+, BH4 , and NH2" ions is essentially an amalgam 



of the structures present in LiNHi and LiBH4. 

Turning now to the dehydrogenation of Li4BN3Hio, we 
note first that although several dehydrogenation products have 
been observed in experimental studieS j"''^'^'^ the precise iden- 
tity and respective proportions of these phases have not yet 
been definitively determinediii For example, Aoki et al. re- 
ported unidentified diffraction peaks after dehydrogenation at 
29 — 25° and 35°. A technique for predicting possible de- 
composition pathways would be of significant value in help- 
ing to understand these reactions. However, decomposition 
pathways and products are difficult to predict a priori. To 
this end, we have scanned through the thermodynamics of a 
large number of candidate dehydrogenation reactions over a 
wide temperature range in order to identify the energetically 
favored products. 

Calculated thermodynamics for seventeen candidate 
Li4BN3Hio dehydrogenation reactions are listed in Table lVIll 
In addition to the room-temperature energetics, there we also 
present data for T — 460 K, which is slightly below the 
melting point of pure Li4BN3Hio. Unlike the Li4BN3Hio 
formation reaction, which can proceed at room temperature^ 
dehydrogenation of Li4BN3Hio has been reported only at 
elevated temperatures. In the case of pure Li4BN3Hio, dehy- 
drogenation occurs above ~520 K, about 60 degrees above 
the melting temperature."' However, a significant reduction in 
dehydrogenation temperature was observed upon the addition 
of a small amount of a PtA'ulcan carbon catalyst.^ In this 
latter case hydrogen release occurred from the solid phase at 
T > 390 K. To more completely assess the dehydrogenation 
thermodynamics over a range of relevant temperatures, in 
Fig.|2]we plot the enthalpies and free energies of the candidate 
reactions for T = 0-1000 K. 

Out of the seventeen reactions considered, three 
reactions — (i), (c), and (f) — emerge as the most-favorable 
reactions in distinct temperature regimes (see Fig. |2] 
bottom panel). For low temperatures reaction (i) is pre- 
ferred, followed at increasing temperatures by reaction (c) 
and reaction (f). More specifically, the favored products 
and their respective temperature ranges of stability are: 

K <T< 300 K: 2LiNH2 + 2LiH H- BN H- 2H2 
300 K <T< 700 K: Li3BN2 + UNHz + 4H2 
700 K <T: Li3BN2 + LiH + iN2 + |H2 

The relatively large entropy of the gas phase products plays 
an important role in determining which reaction is favored at 
a given temperature. Reactions yielding greater quantities of 
gaseous products should be favored with increasing T, and 
this is consistent with the observed trend: reaction (i), 2 mols 

reaction (c), 4 mols reaction (f), 5 mols. (It should be 
noted that at high temperatures — and certainly above the melt- 
ing point of Li4BN3Hio — the harmonic approximation will no 
longer be valid, so we expect some decrease in accuracy with 
increasing T.) 

Of the three favorable dehydrogenation reactions, only re- 
action (c), Li4BN3Hio Li3BN2 + LiNH2 + 4H2, takes place 
within the temperature range where H2-desorption (from un- 
doped Li4BN3Hio) has been experimentally observed (T ^ 
520-630 K). While our prediction of Li3BN2 as a dehy- 
drogenation product is consistent with experimental observa- 
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TABLE VI: Formation energies of Li4BN3Hio [in kJ/(mol Li4BN3Hio)]. AE corresponds to tiie static zero Kelvin DFT energies; Aff^^*"^ 
adds the zero point energies to AE; AH^^ ^"""^ (Eq. [S} adds finite-temperature vibrational (and molecular rotational+translational+pl/) en- 
ergies to A//^"""^; AG^"^™'^ (Eq. includes all of the contributions to Aff^"""""^ and adds vibrational (solids) or tabulated (molecular) 
entropies. 

Rxn.No. Reaction AE AH^^"*" ^^T^aooK ^^t^ISST 

(si) 4U + B + ^-^2 + 5H2 ^ Li4BN3Hio -806.9 -670.2 -708.1 -454.6 

(b) 3 LiNH2 + LiBH4 ^ Li4BN3Hio -11.3 -11.2 -11.8 -10.2 



TABLE VIL Calculated reaction energies [AE in kJ/(mol H2)], enthalpies [Ai7 in kJ/(mol H2)], and free energies [AG in kj/mol products] 
for candidate Li4BN3Hio dehydrogenation reactions. 



Rxn. No. 




Reaction 


AE 




A// 


AH 


AG 


AG 


(c) 


Li4BN3Hio - 


^ Li3BN2 + LiNH2 + 4H2 


28.2 


5.7 


11.2 


12.5 


-88.8 


-161.3 


(d) 


Li4BN3Hio - 


^ Li3BN2 + iLi2NH -1- iNH3-l- 4H2 


40.6 


17.7 


23.4 


24.3 


-61.4 


-145.2 


(e) 


Li4BN3Hio - 


^ Li3BN2 + LiH + NH3 + 3H2 


42.6 


17.5 


23.5 


24.0 


-71.6 


-148.0 


(f) 


Li4BN3Hio - 


Li3BN2 + LiH + + |H2 


50.2 


24.1 


29.7 


30.8 


-38.2 


-131.0 


(g) 


Li4BN3Hio - 


Li3BN2 + Li -H iN2+ 5H2 


61.9 


37.4 


43.5 


44.6 


23.0 


-82.2 


(h) 


Li4BN3Hio - 


U3BN2 +U + NH3 + |H2 


60.4 


37.5 


44.1 


44.7 


-10.5 


-99.2 


(i) 


Li4BN3Hio - 


^ 2LiNH2 + 2LiH -H BN -1- 2H2 


6.4 


-19.4 


-16.1 


-15.4 


-90.1 


-121.2 


Q) 


Li4BN3Hio - 


-* LiNH2 + Li3N -H BN -1- 4H2 


56.9 


35.8 


40.5 


41.3 


34.1 


-34.9 


(k) 


Li4BN3Hio - 


^ 4LiH -1- BN -1- N2 + 3H2 


79.6 


44.1 


48.5 


48.8 


11.1 


-60.7 


(1) 


Li4BN3Hio - 


^ 2Li2NH BN -H 4H2 


45.3 


23.3 


27.9 


28.8 


-14.9 


-83.3 


(m) 


Li4BN3Hio - 


+ 2LiNH2 BN 2Li -1- 3H2 


60.0 


39.5 


45.1 


46.2 


32.3 


-23.6 


(n) 


Li4BN3Hio - 


LiNH2 + Li2NH + LiH -H BN -1- 3H2 


32.3 


9.1 


13.2 


14.1 


-52.5 


-102.3 


(0) 


Li4BN3Hio - 


^ Li2NH -1- 2LiH BN -1- NH3 + 2H2 


55.9 


28.5 


32.7 


32.1 


-35.4 


-88.9 


(P) 


Li4BN3Hio - 


^ LiNH2 + 3LiH BN -1- NH3 + H2 


27.7 


-9.1 


-6.6 


-8.8 


-73.0 


-107.9 


(q) 


Li4BN3Hio - 


^ Li3N -1- LiH -1- BN NH3 + 3H2 


80.9 


57.7 


62.6 


62.4 


51.3 


-21.6 


(r) 


Li4BN3Hio - 


iLi2NH -1- Li3N BN iNH3 + 4H2 


69.3 


47.8 


52.7 


53.2 


61.5 


-18.8 


(s) 


Li4BN3Hio - 


^ Li3BN2 + iLi3N -1- 2NH3 + 4H2 


48.7 


26.0 


31.7 


32.5 


-35.9 


-123.6 



tions, we are not aware of any reports indicating the pres- 
ence of LiNH2. This discrepancy is likely due to the fact that 
most experiments have been performed on off-stoichiometric 
Li4BN3Hi() formed from 2:1 LiNH2:LiBH4 mixtures. Mix- 
tures with this composition are deficient in the Li and N pre- 
sumably needed to yield the additional amide present in our 
calculations. 

The calculated enthalpy for reaction (c) ranges from 11.2 
to 12.8 kJ/(mol H2) for T = 300-500 K, indicating that H2 
release is a weakly endothermic process. For the higher- 
temperature reaction (f), dehydrogenation is likewise pre- 
dicted to be endothermic with a slightly larger enthalpy of 
^31 kJ/(mol H2). Only the low temperature reaction (i) 
with AH w -(16-19) kJ/(mol H2) is found to be exother- 
mic. It would be interesting to test whether the predicted 
temperature-dependent decomposition sequence [(i) (c) 
(f)] could be observed experimentally by holding Li4BN3Hio 
at temperatures within each of the reaction regimes for long 
times. 

Based on the data presented in Fig.|2]and Table lVlH H2 des- 
orption from Li4BN3Hio is thermodynamically favorable at 
essentially all non-zero temperatures: the negative free ener- 
gies of dehydrogenation suggest that Li4BN3Hi() is metastable 
with respect to decomposition via one of the reactions (i), (c), 
or (f). On the other hand, experiments'* on pure Li4BN3Hio 
indicate that desorption does not take place until the tempera- 




400 600 

Temperature (K) 



FIG. 2: (Color online) Calculated enthalpies (AH, top panel) and 
free energies (AG, bottom panel) of candidate Li4BN3Hio dehydro- 
genation reactions as a function of temperature at p = 1 bar. The 
reactions are labeled as in TablelVIll 
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ture is raised above 520 K. Taken together, these factors sug- 
gest that the relatively high temperatures needed in practice 
for Ha-release are a consequence of poor kinetics, not unfa- 
vorable thermodynamics. 

An important distinction between our calculations and the 
experimental dehydrogenation of undoped Li4BN3Hio con- 
cerns the phase of Li4BN3Hi() present during the dehydro- 
genation reaction. As mentioned above, pure Li4BN3Hi() re- 
leases hydrogen from the molten state.'* Similarly, H2 release 
from other hydride phases, such as LiBH4, also occurs from 
the molten state. Due to the computational expense associ- 
ated with obtaining precise energetics for liquids using DFT, 
this study, and previous studies on LiBH4,'*'''*^ have modeled 
the dehydrogenation reaction as a solid-state reaction. Since 
the enthalpy of the liquid phase will be more positive than 
that of the corresponding solid phase, the calculated solid- 
state enthalpy of reaction (c) evaluated at the melting point 
(Tmp), AH'^'^r ~ 12.5 kJ/(mol H2), represents an upper bound 
on the enthalpy for desorption from liquid Li4BN3Hio. It 
will be possible to more precisely estimate the magnitude by 
which /S.H is overestimated once the latent heat of melting 
for Li4BN3Hio has been measured. Until such data become 
available, it nevertheless seems reasonable to conclude that 
the dehydrogenation of Li4BN3Hio is either weakly endother- 
mic, or exothermic. This conclusion is consistent with the 
failed attempts at re-hydriding Li4BN3Hio reported thus far 
in the literature,'* as a small (positive) enthalpy of dehydro- 
genation will result in a similarly small thermodynamic driv- 
ing force for re-hydriding. On the other hand, for Li4BN3Hio 
doped with Pt/Vulcan carbon,^ dehydrogenation begins below 
Tjnp. In this case our approach of modeling solid state re- 
actions more accurately captures the true phase behavior of 
Li4BN3Hio, neglecting the possible formation of C- or Pt- 
containing compounds. 

V. CONCLUSION 

Density functional theory calculations have been employed 
to study the crystal structure and the finite-temperature for- 
mation and dehydrogenation thermodynamics of Li4BN3Hio. 
The calculations have resolved the discrepancies in hydrogen 
bond lengths between the separate structures determined via 
X-ray and neutron diffraction experiments, and suggest that 
the neutron data yields a slightly more accurate description of 
the crystal structure. All three of the reported experimental 



crystal structures relax to essentially the same structure. 

For reaction energetics, our calculations indicate that 
Li4BN3Hio formation is exothermic, and that the dehydro- 
genation enthalpy of solid-state Li4BN3Hio is temperature- 
dependent: H2 release is exothermic at low temperatures, and 
weakly endothermic at the higher temperatures probed by re- 
cent experiments. A non-zero latent heat for Li4BN3Hi() melt- 
ing will reduce the dehydrogenation enthalpy from the liquid 
state (relative to the solid state), but the size of this effect 
is still to be determined. To help clarify the identity of the 
phases produced via dehydrogenation we have performed a 
computational search over 17 candidate dehydrogenation re- 
actions, and identified three reactions having favorable ther- 
modynamics spanning the temperature range T — 0-1000 K. 
All three of these reactions exhibit a decrease in free energy, 
suggesting that Li4BN3Hi() is a metastable phase. For tem- 
peratures where H2-desorption has been experimentally ob- 
served, the thermodynamically-favoredreaction is Li4BN3Hio 
Li3BN2 + LiNH2 + 4H2, with an enthalpy of 11.2-12.8 
kJ/(mol H2). The relatively small dehydriding enthalpies are 
consistent with the failed re-hydriding attempts reported in the 
literature, and suggest that hydrogen release from Li4BN3Hio 
is a kinetically — rather than thermodynamically — hindered 
process. 

Note added in proof — During the review of this manuscript 
two new studies of Li4BN3Hio were brought to our 
attention. In the first study Noritake and co-workers^" per- 
formed a crystal structure analysis of Li4BN3Hi() using syn- 
chrotron X-ray diffraction. While they found external struc- 
tural parameters (bcc, space group /2i3, a — 10.67 A) con- 
sistent with the two reports cited previouslyr^i^ the reported 
N-H bond lengths were more consistent with the neutron data 
of Ref. 12 (and with our calculated N-H distances) than the 
X-ray data of Ref. [ll]. In the other recent study, Herbst and 
Hector"*^ examined the electronic structure and energetics of 
of Li4BN3Hi() using DFT. Our results appear to be consistent 
with their findings. 
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modynamic properties — for example, ZPEs calculated with both 
PAW sets differ by 0.7 kJ/mol f.u. at most. It therefore seems 
likely that an approach where all vibrational properties (solids and 
molecules) were evaluated using the softer potentials would also 
yield reasonably precise results. We thus conclude that our ap- 
proach of mixing potentials gives results comparable in precision 
to those obtained with the hard PAW set alone, but without the 
associated computational expense. 

The discrepancy in calculated vs. experimental enthalpy for 
LiiNH has been noted in Ref. and it has been suggested that 
the experimental measurement be revisited in light of the good 
agreement for the other Li-containing compounds: LiH, LiNH2, 
and LiBH4. To our knowledge the formation enthalpy of Li2NH 
reported in Ref. 46 has not been independently confirmed. 
Due to the computational difficulty associated with treating a 
non-stoichiometric phase, our calculations are restricted to the 
stoichiometric composition "Li4BN3Hio," whereas most exper- 
iments have focused on 2:1 LiNH2:LiBH4 off-stoichiometric 
"Li3BN2H8" mixtures, which have been found to minimize^ NH3 
release. 



